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We examine the spectrum of superparticles obtained from the minimal 
50(10) grand unified model, where it is assumed the gauge symmetry break- 
On. ing yields the Minimal Super symmetric Standard Model (MSSM) as the ef- 
r~| ! fective theory at Mgut ~ 2 x 10 16 GeV. In this model, unification of Yukawa 
couplings implies a value of tan /3 ~ 45 — 55. At such high values of tan j3, 
assuming universality of scalar masses, the usual mechanism of radiative elec- 
troweak symmetry breaking breaks down. We show that a set of weak scale 
sparticle masses consistent with radiative electroweak symmetry breaking can 
be generated by imposing non-universal GUT scale scalar masses consistent 
with universality within 50(10) plus extra .D-term contributions associated 
with the reduction in rank of the gauge symmetry group when 5O(10) sponta- 
neously breaks to SU (3) x SU(2) x [7(1) . We comment upon the consequences 
of the sparticle mass spectrum for collider searches for super symmetry. One 
implication of 5O(10) unification is that the light bottom squark can be by 
far the lightest of the squarks. This motivates a dedicated search for bottom 
squark pair production at pp and e + e~ colliders. 

PACS numbers: 14.80.Ly, 13.85. Qk, 11.30.Pb 
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Unification of the Standard Model (SM) of strong, weak and electromagnetic interactions 
within a single Lie group such as SU(5) or 50(10) has a long history and many attractive 
features [|TJ. £77(5) is the smallest grand unifying group, and predicts the quantization of 
electric charge, the unification of gauge couplings and the unification of bottom and tau 
Yukawa couplings at scales of Q = Mqut — 10 15 GeV [0. The 50(10) theory incorporates 
all the matter fields of the SM into the 16-dimensional spinor representation, x/jxq, of 50(10) 
||. In minimal 5O(10), not only the gauge couplings but all the Yukawa couplings (within 
a generation) are unified at Q = Mqut- If the right-handed neutrino field present in ip ie 
acquires a large Majorana mass, it decouples from the theory, and a small neutrino mass is 
induced via the see-saw mechanism [|J. 

The supersymmetric version of this model, with supersymmetry (SUSY) softly broken 
at a scale ^ 1 TeV naturally stabilizes the hierarchy between the weak scale and the grand- 
unification scale. Supersymmetry also raises the unification scale to Mqut — 2 x 10 16 GeV, 
which helps reduce the rate for proton decay to below the level of experimental bounds. 
In addition, the introduction of supersymmetry with soft SUSY breaking (SSB) masses of 
order the weak scale allows for the near unification of gauge coupling constants |J. In 
supergravity-based models, it is usually assumed that all scalar masses receive a common 
mass niQ = mu = = m L = the = rriH u = m H d = m o at M GUT , while all gauginos 
receive a common mass m 1( / 2 and all trilinear SSB terms unify to A . The SSB masses 
and couplings are then evolved via renormalization group equations (RGEs) from Mqut to 
Q ~ M wea k. The m\ term is driven to negative values, which results in radiative breaking 
of electroweak symmetry, provided the top quark mass is large (e.g. 175 GeV). 

In addition to the matter superfield ip±6, the minimal 5O(10) model includes a 10 di- 
mensional Higgs superfield </> 10 that decomposes into a 5 + 5 representation of £77(5), and 
includes the two Higgs superfields (H u and H d ) of the Minimal Supersymmetric Standard 
Model (MSSM). The superpotential includes the term 

W 3 + • • ■ 

responsible for quark and lepton masses, with A the single Yukawa coupling in the low 
energy theory. The dots represent terms including for instance higher dimensional Higgs 
representations and interactions responsible for the breaking of £O(10). 

The mass spectrum of SUSY particles in minimal supersymmetric £O(10) constrained 
by radiative electroweak symmetry breaking has been studied previously in a number of 
papers [|6|-|T6|j. Unification of bottom, tau and top Yukawa couplings was found to occur 
at very large values of the parameter tan (3 ~ 50 — 60, and specific spectra were generated 
for values of m t ~ 190 GeV 0. Assuming universality of soft SUSY breaking masses at 
Mqut, it was found ^,TD] that Yukawa unification consistent with radiative electroweak 
symmetry breaking could also occur for m t < 170 GeV as long as m^j ^ 300 GeV. This 
generally leads to sparticle masses far beyond the reach of the CERN LEP2 or Fermilab 
Tevatron pp colliders. For values of mt — 175 GeV, solutions including radiative electroweak 



breaking were very difficult to achieve. In Ref. | lq| , the SUSY particle mass spectrum was 
investigated with non-universal SSB masses. Various solutions were found, but the non- 
universality in general broke the £O(10) symmetry. In Ref. [Tj|, it was argued that £O(10) 
D-term contributions to scalar masses had the correct form to allow for successful radiative 
electroweak symmetry breaking and the computation of weak scale SUSY particle masses. 
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In this report, we explicitly calculate the sparticle mass spectrum for SO (10) SUSY 
GUT models, taking the pole mass m t = 175 GeV. We make the following assumptions. We 
assume the structure of minimal SUSY SO (10) above the scale Q = M GUT . We assume that 
SUSY 50(10) directly breaks to the MSSM at M GUT . Accordingly, there exist independent 
masses mie and mio for the matter and Higgs scalar fields. In the breakdown of SO(10) to 
SU(3)c x SU(2)l x U(l)y, additional D-term contributions (parametrized by M|, which 
can be either positive or negative) to the SSB scalar masses arise JT7[ : 

m Q= m E = m U = m 16 + M D 

m 2 D = m 2 L = m 2 6 — 3M D 
m k, d = m ioT2M 2 D . 



Thus, the model is characterized by the following free parameters: 

mie, m 10 , Mf,, mi/2, A , sign(n). 

The value of tan will be restricted by the requirement of Yukawa coupling unification, and 
so is not a free parameter. 

Our procedure is as follows. We generate random samples of model parameters 

< mie < 1500 GeV, 

< mio < 1500 GeV, 

< m V 2 < 500 GeV, 
-500 2 <M 2 D < +500 2 GeV 2 , 

45 < tan (3 < 55, 
-3000 <A < 3000 GeV and 

[A > or /i < 0. 

We then calculate the non-universal scalar masses according to formulae given above, and 
enter the parameters into the computer program ISASUGRA. ISASUGRA is a part of the 
ISAJET package |18[] which calculates an iterative solution to the 26 coupled RGEs of the 
MSSM. 

To calculate the values of the Yukawa couplings at scale Q = Mz, we begin with the pole 
masses 77i& = 4.9 GeV and m T = 1.784 GeV. We calculate the corresponding running masses 
in the MS scheme, and evolve and m T up to M z using 2-loop SM RGEs. At Q = M z , we 
include the SUSY loop corrections to and m T using the approximate formulae of Pierce 



et al. |i~9|l . A similar procedure is used to calculate the top quark Yukawa coupling at scale 
Q = m t . 

Starting with the three gauge couplings and t, b and r Yukawa couplings of the MSSM at 
scale Q = Mz (or m t ), ISASUGRA evolves the various couplings up in energy until the scale 
where g± = g2, which is identified as Mqut, is reached. The GUT scale boundary conditions 
are imposed, and the full set of 26 RGE's for gauge couplings, Yukawa couplings and relevant 
SSB masses are evolved down to Q ~ M wea k, where the renormalization group improved one- 
loop effective potential is minimized at an optimized scale choice Q = ■yJfnf L m^ R and radiative 
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electroweak symmetry breaking is imposed. Using the new spectrum, the full set of SSB 
masses and couplings are evolved back up to M GUT including weak scale sparticle threshold 
corrections to gauge couplings. The process is repeated iteratively until a stable solution 
within tolerances is achieved. We accept only solutions for which the Yukawa couplings X t , 
Xb and A r unify to within 5%. This constraint effectively fixes the value of tan/3 typically 
to ~ 48. Yukawa unified solutions are found only for values of /J, < 0. We also require 
the lightest SUSY particle to be the lightest neutralino, and that electroweak symmetry is 
successfully broken radiatively. 

We show in Fig. |I]the regions of model parameter space for which a SUSY mass spectrum 
can be calculated consistent with the above constraints. In Fig. [L]a, we show the plane of 
mio vs. rriie. Each dot represents a point for which a solution was obtained. Points denoted 
by a cross are valid solutions, but with sparticle or Higgs masses below existing limits from 



LEP2. We require m fl > 73 GeV, m~± > 95 GeV and m h > 85.2 GeV f2(|. From the 
distribution of points, we see that regions of model parameter space with m 16 < m 10 are 
preferred, although for very large values of mi 6 , a few solutions are obtained for m w > mi 6 . 
In Fig. |1]&, we plot the vs. ttliq parameter plane. In this frame, actually stands for 
sign(Ml)) x J\M^\. No solutions were obtained for M|> < 0, and in fact no solutions were 



obtained for M|, = 0: this illustrates that non-zero D-term contributions to scalar masses are 
crucial for a valid sparticle mass spectrum in minimal 5*0(10). The requirement of positive 
definite D-term contributions to scalar masses will leave, as we shall see, a distinctive imprint 
on the SUSY particle mass spectrum [[17]]. From the mi/2 vs. m 16 plane in Fig. |l]c, it can be 
seen that m\§ is typically larger than m^; otherwise f\ becomes the lightest SUSY particle, 
in violation of cosmological limits on charged relic particles. Finally, in Fig. we show 
the range of tan/3 values for which solutions were generated versus the parameter rrii/2- We 
see that 46 < tan (3 < 52, with the slightly higher values of tan j3 being preferred when m^j 
is large. The bounds on tan j3 are weakened if r — b — t Yukawa unification is relaxed to 
more than 5%. 

In Fig. we show the range of selected sparticle and Higgs boson masses that are 
generated within minimal SO(10) with Yukawa coupling unification. In frame a), we see 
that the light Higgs boson h has mass generally bounded by rrih < 125 GeV. This range of 
light Higgs boson masses may well be accessible to Fermilab Tevatron Higgs boson searches 
21] . Values of rrih ^110 GeV are associated with cases where rriA becomes comparable to 



or smaller than Mz- In frame 6), we plot solutions in the /i vs. M 2 plane, where M 2 is the 
SU(2) gaugino mass. Many solutions with < M 2 exist, which generally implies that the 
lighter charginos and neutralinos have substantial higgsino components. The solutions with 
large \fj,\ and M 2 ~ 100 GeV all correspond to values of rriiQ > 1300 GeV. In frame c), the 
bottom squark mass mi is plotted versus rna R . We see that although m ifl can be only as 
light as ~ 700 GeV, the b\ mass can be as low as ~ 150 GeV. The bottom squark (mainly 
6r) is generically much lighter than other squarks, because of the D-term contribution to 
mn at Q = Mqut as well as b- Yukawa coupling effects which are significant for large values 
of tan f3. Finally, in frame d), we show the lightest tau slepton mass versus the light chargino 
mass. In S'O(IO), the stau is the lightest of the sleptons, but as can be seen, solutions with 
< 200 GeV are very difficult to generate, and almost always, > rn~±, so that two 

body decays such as x 2 — > f±r or xf ~^ T\V almost never occur. A feature of minimal 
SO (10) is that the light stau may contains a large left stau component, whereas in models 
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with universality, the light stau is dominantly a right slepton. This could have an impact 
on the efficiency of detecting daughter tau leptons via their hadronic decay. 

In Table I, we show sample weak scale sparticle and Higgs boson masses for five SO (10) 
solutions with unified Yukawa couplings. It is possible to find solutions with sparticle 
masses potentially accessible to both LEP2 and Fermilab Tevatron searches, in contrast 
to previous studies assuming universality of scalar masses at Mqvt- For case 1, we take 
(mi 6 , mio, M D , mi/a, A ) = (405.8, 680.3, 96.8, 427.2, 596) GeV. This solution requires 
tan/5 = 51.3 to unify the Yukawa couplings. The evolution of gauge and Yukawa couplings 
for this case is shown in Fig. |3|a. In our program, we do not require the SU(3) gauge 
coupling to exactly unify with the SU(2) and U(l) gauge couplings, but rather attribute the 
near miss to unknown high scale physics. The Yukawa couplings diverge from their unifica- 
tion point and evolve to M wea k, with a kink in the curves coming from weak scale threshold 
effects. In Fig. we show the evolution of SSB Higgs boson masses and third generation 
SSB masses. We actually plot sign(m 2 H ) x \J\m 2 H \. In this case, both Higgs squared masses 
evolve to negative values, signaling the onset of radiative electroweak symmetry breaking. 
The GUT scale non-universality due to D-term contributions is evident. It usually results 
in left SSB slepton masses being close to or lighter than right slepton masses, and right 
sbottom masses lighter than the other squark masses. 

The final weak scale sparticle masses are listed in Table I. For case 1, none of the 
sparticle or Higgs bosons are accessible to LEP2, while one or more of the Higgs bosons may 
be accessible to the Fermilab Tevatron running at maximal luminosity. An e + e~ collider 
operating at yfs = 500 GeV would find not only the various MSSM Higgs bosons, but also 
charginos and neutralinos (with substantial higgsino components) and light tau sleptons (in 
this case f\ ~ tr). 

The second case study point shown has very large values of m\§ and mio, leading in 
general to a spectrum with very heavy scalars. The exception in this case is that the b\ 
mass is only 140 GeV, and is directly accessible to Fermilab Tevatron collider searches, 
according to Ref. [p3| . 

In this case, b% — > bxi with a branching fraction of ~ 80%, so that bibi production would 
be visible in bb+ fix events. The gluino is also relatively light and decays via g — > bb\. 
it would be interesting to examine whether the improved b tagging at Tevatron upgrades 
would allow its detection in the multi-f> plus fix channel. The light Higgs boson might 
be accessible to high luminosity Tevatron experiments, but the charginos and neutralinos 
would be difficult to see via the trilepton channel since X2 dominantly decays to bbxi and 
the X2 ~* e ^Xi branching fraction is only 0.8%. A Higgs signal of bb£+ $t events from Wh 
production would contain substantial contamination from X1X2 events, which give rise to 
the same event topology. Contrary to models with universality, the left selectron (smuon) is 
significantly lighter than the right selectron (smuon). This distinctive feature of the SO (10) 
model would be difficult to discern as the sleptons are very heavy. In the squark sector, the 
cIr and sr are the lightest of the first two generations of squarks, owing to the D-terms. 

In case 3, again the bulk of the scalars are quite heavy, and well beyond the reach of LEP2 
or the Tevatron. Again the exception is the light bottom squark. In this case, however, b\ 
decays with a 24% (8%) branching fraction to bx\ so the event signatures will be more 

complicated. Since X2 an d X3 decay with a large rate to bbxi, some of the b\bi events will 
contain final states with up to six 6-jets plus If clean trilepton signatures are detected 
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2~4| , they will contain a mixture of events from both X1X2 an d X1X3 production. 

In case 4, all strongly interacting sparticles including the bottom squark are quite heavy 
and accessible only at the LHC However, the various sleptons and sneutrinos are within 
reach of an e + e~ linear collider operating at y/s ~ 1000 GeV. In this case, a very mixed fi 
whose composition may be measurable at a Linear Collider may serve to distinguish this 
framework from models with universal soft masses. Moreover, z/ e ^ and the are measurably 
lighter [273] than e.R, again in contrast with expectations in models with universality. Note 



also that < M 2 , so that the light charginos and neutralinos have substantial Higgsino 
components, and further that there is only a small mass gap between m~o or m~± and m~o, 
so that -ino decay products will be soft. 

Finally, in case 5, again the light charginos and neutralinos are higgsino- like, and will be 
challenging to detect at the Fermilab Tevatron collider. This spectrum is characterized by a 
very light Higgs boson spectrum, and in fact 32% of top quark decays are to charged Higgs 
bosons. Indeed, both the light and pseudoscalar Higgs boson are at the edge of detectability 
at the LEP2 collider. 

We have demonstrated that the inclusion of D-terms can lead to radiative electroweak 
symmetry breaking even in models with Yukawa coupling unification. As shown in Fig. |l]&, 
we are unable to find corresponding solutions for models with scalar mass universality for 
the ranges of parameters studied here. We have not attempted to do an analysis of the 
phenomenological implications of the model. In a follow-up report, we will present results of 
calculations for the neutralino relic density, b — *> 57 decay rate, direct dark matter detection 
rate, and prospects for collider searches [^7]]. Parts of the parameter space as well as some 
of the case studies may well be excluded by experimental constraints. For instance, our 
preliminary results indicate that the predicted value for the decay b — > 57 exceeds the 
experimental upper limit by a factor ~ 2 — 4 if m^j ~ 200 — 500 GeV. This is well-known 



to be a problem common to models with large tan (3 and fi < |28| , the region of parameter 
space where Yukawa couplings unify. 

Despite this phenomenological problem, we find it encouraging that it is possible to 
construct a calculable framework with gauge and Yukawa coupling unification. We can 
imagine other physics that may make it possible to circumvent experimental limits such 
as those from b — > 57. For example, it has been pointed out |29| that if the right-handed 
neutrino mass is significantly below the GUT scale and if mf = 2m^ 6 , third generation 
scalars would be radiatively driven to much lower masses than other matter scalars, and 
further, that when D-terms are included, radiative electroweak symmetry breaking is still 
possible |30| . In this case, since the degeneracy between squarks is badly broken, it is possible 



that gluino-mediated contributions to b — > 57 (which are generally thought to be small) 
may be significant. Whether these are large enough (and of the correct sign) to cancel the 
chargino-mediated amplitudes remains to be investigated. Alternatively, one might imagine 
that the usual computation of b — > 57 amplitudes may be altered by CP violating phases 
between various chargino amplitudes , or by large Yukawa coupling radiative corrections 



Summary: We have shown that explicit evaluation of sparticle mass spectra is possible 
in the minimal SUSY SO(10) model with Yukawa coupling unification and radiative elec- 
troweak symmetry breaking, by including non-universal SSB masses at Q = Mqut which are 
in accord with SO(10) breaking to the gauge group of the MSSM. The resulting spectra re- 
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fleet the influence of the D-term contributions to scalar masses. Characteristic features of the 
model can include a light sbottom, a f\ which is mainly ?l, < and lighter charginos 
and neutralinos with substantial (sometimes even dominant) Higgsino components. 
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TABLES 



TABLE I. Weak scale sparticle masses and parameters (GeV) for five 50(10) case studies. 



parameter 


case 1 


case 2 


case 3 


case 4 


case 5 


mi 6 


405.8 


1240.0 


1022.0 


414.8 


629.8 


m w 


680.3 


1414.0 


1315.0 


735.7 


836.2 


M D 


96.8 


410.6 


329.8 


171.9 


135.6 


m 1/2 


427.2 


136.5 


232.0 


449.1 


348.8 


A 


596.0 


-1100.0 


-1350.0 


57b. 7 


-186.5 


tan p 


Dl.O 


4/.U 


AQC 

40.0 


ol.o 


£9 1 


m g 


lUZl.o 


/i no n 
4U9.9 


631.5 


1U69.3 


8o4.o 


mu L 


9o3. I 


133 / .8 


11/8.5 


1U33. / 


y /4.4 


m 3 


09^ /I 


1U0 1 .0 


Q7n 1 




Q1 n s 
y lu.o 




718 4 


737 9 


512 3 


754 5 


618 7 


mr 

01 


735.6 


140.6 


187.1 


721.7 


636.8 


m,- 


478.7 


1012.8 


857.8 


428.3 


634.6 


m« 
£r 


452.9 


1321.1 


1088.9 


489.6 


662.5 


472.0 


1009.7 


854.1 


420.7 


629.5 




233.2 


790.1 


623.6 


272.6 


427.8 




386.0 


787.4 


619.5 


314.6 


519.1 


m x% 
m x\ 


159.5 


110.5 


122.9 


177.5 


106.3 


166.7 


110.3 


131.6 


195.1 


126.1 


129.0 


56.8 


84.0 


152.3 


87.5 


m h 


113.7 


115.5 


118.8 


116.4 


93.7 


m A 


115.8 


645.0 


479.9 


277.9 


93.9 


m H + 


152.2 


652.3 


490.2 


295.1 


137.1 




-157.2 


-329.8 


-150.5 


-185.5 


-113.9 


(ti\tl) 


0.14 


0.99 


0.99 


0.47 


0.11 
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FIG. 1. Plots of regions of parameter space where valid solutions to minimal SUSY SO(10) are 
obtained, consistent with Yukawa coupling unification to 5%, and radiative electroweak symmetry 
breaking. 
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FIG. 2. The range of selected sparticle masses that are generated in minimal SUSY S0(W) 
models with Yukawa coupling unification and radiative electroweak symmetry breaking. 
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FIG. 3. For case 1 in Table I, we show a) the running of both gauge and Yukawa couplings 
between Q = Mqut and Q = M wea k. In b), we show the running of SSB Higgs masses (dashed 
curves) and third generation SSB masses (solid curves). 
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